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Figure  11.  Tidal  rest  currents  in  Strait  of  Gibraltar 
(HN  model). 

Figure  12.  Measured  rest  currents  in  Strait  of  Gibraltar 
(after  Lacombe,  1961). 

Figure  13.  Measured  rest  currents  at  different  depths  in 
Strait  of  Gibraltar  (locations,  see  Figure  12) 
(after  Lacombe ,  1961). 

Figure  14.  Transport  of  water  through  Strait  of  Gibraltar 
from  surface  to  100m  (HN  model). 
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Abstract 

The  Hydrodynamical  Numerical  (HN)  model  of  Walter  Hansen 
has  been  used  to  compute  tides  and  currents  in  the  Strait  of 
Gibraltar.  The  model  has  two  open  boundaries  at  which  the  tides 
were  prescribed  at  each  time  step,  using  four  tidal  constituents. 
The  grid  size  was  two  nautical  miles,  and  the  time  step  (from 
Courant  criterion)  was  24  seconds.  Equilibrium  was  established 
after  20  hours  of  real-time  computation.  The  main  results  from 
these  model  runs  were:  (11  At  mean  tides,  the  inflowing  currents 
are  stronger  and  last  somewhat  longer  than  the  outflowing  currents. 

(2)  There  is  a  net  inflow  into  the  Mediterranean  Sea  caused  by 
tidal  currents  alone  (excluding  wind  and  sea  level  difference). 

(3)  The  tides  computed  with  the  HN  model  were  in  good  agreement 
with  those  obtained  from  the  harmonic  method  at  points  where 
harmonic  constants  were  available.  There  was  also  a  good  agree¬ 
ment  between  computed  currents  and  the  few  available  current 
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measurements.  (41  The  HN  model  gave  a  net  transport  of  1  »48  km 

.< 

per  hour  into  the  Mediterranean  in  the  upper  100  m.  This  value  is 
somewhat  higher  than  previous  estimates;  however,  heat  budget 
considerations  require  a  higher  water  exchange  through  the  Strait 


of  Gibraltar  than  do  previous  estimates. 
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1.  INTRODUCTION  (Purpose  of  the  paper). 

The  Hydrodynamical  Numerical  models  of  Professor  Walter  Hansen, 
University  of  Hamburg,  have  been  tested  on  a  number  of  shallow  semi- 
closed  seas.  The  model  closely  reproduces  observed  tides,  storm 
surges  and  currents  (Hansen,  1966;  SUndermann,  1966;  Laevastu  and 
Stevens,  1969).  It  was  desirable  to  test  this  model  on  the  Strait  of 
Gibraltar,  which  is  a  relatively  deep  area  and  has  two  open  boundaries 
with  distinctly  different  tides.  It  is  of  interest  to  know  whether  the 
latter  aspect  is  properly  reproduced  by  HN  models.  Furthermore, 
the  Strait  of  Gibraltar  is  an  ideal  area  for  eventual  testing  of  a  two- 
layer  HN  model.  This  report  presents  the  results  and  some  verifica¬ 
tions  of  single-layer  HN  model  tests  for  the  Strait  of  Gibraltar. 

2.  HYDRODYNAMICAL  FORMULAS  AND  W.  HANSEN'S  FINITE 

DIFFERENCE  METHOD. 

The  hydrodynamical  numerical  (HN)  method  for  computation 
of  currents  and  sea  level  changes  was  proposed  in  its  present  form 
by  Hansen  in  1956 . 

The  following  basic  equations  are  used  in  the  single-layer 
model: 


(and  )  are  usually  expressed  as: 
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yk)  =  v/ w2  +  w2 

'  ^  x  V  x  y 

The  bottom  stress  (friction)  term  in  formulas  1  and  2  is: 

JQo)  _  _JL_  \f  2  ,  2  .  _r_ 

•  h  u  r  u  v  '  h 

The  following  symbols  were  used  in  the  formulas  above: 
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space  coordinates 
time 

components  of  velocity 

total  depth  (H  =  h  +^) 

surface  elevation 

components  of  external  forces 

components  of  wind  stress  (^  =  3.5x10  ) 

acceleration  of  gravity 

Coriolis  parameter 
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friction  coefficient  (3  x  10  )  (bottom  stress) 
coefficient  of  horizontal  eddy  viscosity 
Laplace  operator 

coefficient  of  friction  (drag  coefficient) 
wind  speed  components 
bottom  stress 


(4) 


(5) 


Analytical  solution (s)  to  formulas  1  to  3  are  of  little  value,  as 


exact  solutions  are  possible  only  for  basins  of  regular  shape,  simple 


depth  and  simple  wind  distribution.  However,  the  formulas  can  be 


solved  on  computers  using  "step-by-step"  finite-difference  methods. 
^t+^(n,m)  ^  (n,m)  -  y-  £  (n,m)  U^n.ml-H^  (n.m-DU^n.m-l)  + 

+  H*  (n-JjlV*  (n-1  ,m)-H*  (n,m)  Vt  (n,m)  ^  (6) 

Ut+2^*(n,m)  =  ^1  -  /H*  (n , m jj U  1  (n , m) 2  +V  *  (n , m )J"  U  (n,m)  + 

+2T*fV  l(n,m)  -  t+^(n,m+ll  -%  t+^(n,m)^  + 

+2  <Xt+2^(n,m)  (7) 

V1  (n,m)  =  2^~(n , t {n, m)2  +  U  V1  (n,m)  - 

-27*fU  t  (n,m)  -  ££  t+  (n,m)-^  t+^(n+l,m)^ + 

+27'Yt+2T\n,m)  (8) 

The  “averaged"  velocity  and  water  elevation  (sea  level)  components  are: 


U  *  (n ,  m)  =#CUt(n,m)  + 


^U^n-l.m)  +Ut  (n+l,m)+Ut  (n^m+D+U^  (n,m-l)^ 
V^n^m)  and  ^t(n,m)  are  analogous. 


(9) 


(The  factor  can  be  interpreted  as  "horizontal  viscosity  parameter." 

Its  normal  value  is  0.99  (see  further  Chapter  4.1)) 

U*t(n,m)=~-  £ Ut(n,m-l)+Ut(n+l,m-l)+Ut(n,m)+Ut(n+l,m) ^  (10) 

*t  *t 

V  (n,m)  is  analogous  to  the  U  (n,m)  above. 
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The  time  step  is  2  T* .  The  total  depth  (H^/H^)  is  computed  as  ' 

(n ,  m)  =  hu  (n ,  m)  +  -y  t+^(n ,  m)  t+T*(n ,  m+  1)J  (11)  r 

The  effects  of  wind  (external  force)  are  computed  with  the  following 
formula: 


The  effects  of  wind  (external  force)  are  computed  with  the  following 
formula: 


A  number  of  slightly  modified  finite  difference  schemes  are 
possible  for  solving  the  hydrodynamical  equations.  Some  of  these 
schemes  are  being  tested  at  FNWC. 
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3 .  HYDRODYNAMICAL  NUMERICAL  (HN)  MODEL  OF 

STRAIT  OF  GIBRALTAR. 

The  computational  grid  for  the  Strait  of  Gibraltar  has  a  two 
nautical  mile  mesh  length.  This  grid  size  and  the  maximum  depth 
in  the  area  require  (according  to  Courant,  Friedrich,  Lewy  criterion) 
a  time  step  of  24  seconds  in  order  that  the  computations  remain 
stable.  Depths  at  u  and  v  grid  points  were  obtained  from  navigation 
charts.  The  friction  coefficient  was  0.03.  Tides  with  four  tidal 
constituents  (M^ ,  S^ ,  N^  and  )  were  introduced  at  each  time  step 
at  both  ends.  The  tidal  input  was  constant  across  the  openings. 

It  has  been  found  at  FNWC  that  the  tidal  input  at  the  openings  can 
be  made  an  inverse  function  of  depth.  However,  if  this  is  not  done, 
the  program  will  adjust  the  z  and  u,  v  values  a  few  gridpoints  inwards. 
Equilibrium  in  this  particular  program  was  established  in  about  20 
hours  of  real  time.  The  program  was  run  100  hours.  Mean  tides 
(tides  between  spring  and  neap)  and  the  hourly  outputs  were  taken 
from  the  last  25  hours  of  computation.  The  program  was  also  run 
with  wind  inputs,  but  essentially  calm  wind  results  are  reported  in 
this  paper. 

4.  VERIFICATION  OF  SEA  LEVEL  COMPUTATIONS. 

The  tides  at  a  number  of  locations  for  which  tid^l  harmonic 
constants  were  available  were  computed  with  the  well-known  harmonic 
method.  The  same  time  span  as  the  output  from  the  HN  method  was 
used  as  well  as  the  same  four  tidal  constituents.  The  sea  level 
changes  computed  with  the  HN  method  were  taken  from  a  grid  point 
closest  to  the  original  location  of  the  corresponding  tide  gauge 
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constituents  in  the  harmonic  computations.  Some  comparisons 
are  shown  in  Figures  1  to  4. 

At  Tangier  (Figure  1) ,  there  is  virtually  no  time  lag  between 
the  HN  and  harmonic  models.  However,  the  HN  model  gives  about 
10%  higher  amplitude,  the  difference  occurring  mainly  in  low  water 
prediction.  At  Tarifa  (Figure  2),  the  tides  with  the  HN  method 
precede  about  10  minutes  the  tides  from  the  harmonic  method  and 
the  HN  tides  are  about  8%  higher.  At  Gibraltar  (Figure  3) ,  and  at 
Ceuta  (Figure  4) ,  the  tidal  heights  of  both  methods  are  in  relatively 
good  agreement;  however,  the  HN  tides  precede  the  "harmonic"  tides 
by  about  30  minutes.  It  should  be  noted  that  the  HN  tides  are  taken 
from  grid  points  which  are  as  much  as  two  miles  away  from  the  tide 
gauges  and  usually  over  deeper  water. 


5.  VERIFICATION  OF  COMPUTED  CURRENTS. 

In  addition  to  producing  tidal  forecasts  in  areas  where 
harmonic  constants  are  not  available,  the  HN  model  has  the  advan¬ 
tages  of  handling  superimposed  meteorological  tides  (storm  surges) 
and  of  predicting  tidal  currents.  The  latter  were  not  obtainable 
using  "harmonic"  methods  unless  current  recordings  were  available 
over  relatively  long  intervals. 

Examples  of  synoptic  tidal  currents  in  the  Strait  of  Gibraltar 
are  given  in  Figures  5  and  6.  The  direct  verification  of  synoptic 
currents  is  difficult  due  to  the  scarcity  of  current  measurements. 

Some  quasi-synoptic  current  charts  have  been  prepared  subjectively 
using  a  variety  of  data,  including  ships'  logbooks.  Gaps  have  often 
been  filled  in  using  some  degree  of  "artistic  license."  Two  subjective 
tidal  current  charts  for  the  Strait  of  Gibraltar  are  given  in  Figures  7 
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and  8  (from  NOO  Publ.  700),  corresponding  to  the  times  and  currents 
in  Figures  5  and  6.  Although  detailed  comparison  of  the  two  sets 
of  charts  has  no  direct  use,  it  can  be  noted  that  agreement  is  good. 

A  more  detailed  comparison  of  measured  and  computed  tidal 
currents  at  a  given  location  in  the  western  entrance  to  the  Strait  of 
Gibraltar  is  given  in  Figure  9.  Although  the  exact  stage  of  tides 
(mean,  spring  or  neap)  during  the  measurement  is  not  known,  the 
relation  between  the  observed  and  computed  tides  is  a  good  one. 

The  tidal  currents,  computed  with  the  HN  method  are  about  20%  stronger. 

It  should  be  noticed  that  the  tidal  currents  are  considerably 
weaker  in  the  Bay  of  Cadiz  than  at  the  entrance  to  and  inside  the 
Strait  of  Gibraltar.  An  example  of  currents  in  the  Bay  of  Cadiz 
computed  with  another  HN  model  is  given  in  Figure  10. 

6 .  TRANSPORT  OF  WATER  THROUGH  THE  STRAIT  OF  GIBRALTAR. 

The  general  circulation  in  the  Strait  of  Gibraltar  (inflow  of 
Atlantic  water  at  the  surface,  outflow  of  Mediterranean  water  along 
the  bottom)  has  been  discussed  in  a  number  of  publications  and 
various  quantitative  estimates  of  the  water  exchange  have  been  given 
(e.g.  Wtlst  1952,  1959;  Lacombe  1961,  and  others).  The  possible 
influence  of  different  forces  in  effecting  the  seasonal  changes  of 
the  water  exchange,  such  as  atmospheric  pressure  difference,  winds, 
difference  in  evaporation  in  the  Mediterranean,  etc.  ,  have  also  been 
discussed.  However,  as  far  as  can  be  ascertained  by  the  authors, 
the  net  (or  rest)  tidal  circulation  in  the  water  exchange  in  the  Strait 
of  Gibraltar  has  not  been  pointed  out  earlier. 


Due  to  the  variation  of  amplitude  and  phase  of  the  tide  over 
some  distance,  there  will  be  a  net  transport  of  tidal  currents.  This 
net  circulation  over  a  mean  tidal  cycle  (24h  50  min)  in  the  Strait 
of  Gibraltar  is  shown  in  Figure  11.  As  seen  from  this  figure,  the 
tides  alone  might  effect  the  net  inflow  to  the  Mediterranean.  This 
is  due  to  the  amplitude  difference  at  opposite  ends  of  the  straits. 

The  best  measurements  of  currents  in  the  Strait  of  Gibraltar 
have  been  made  by  Lacombe  (1961).  The  net  currents,  computed 
from  these  measurements,  are  shown  in  Figures  12  and  13.  As  seen, 
the  main  inflow  occurs  above  100m  depth. 

The  transport  of  water  through  the  straits  from  0  to  100m  during 
a  tidal  cycle  was  extracted  from  the  HN  model  and  is  presented  in 

3 

Figure  14.  The  net  inflow  is  1.43  km  /hour.  This  value  is  about 
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half  of  that  obtained  by  Lacombe  from  direct  measurements^ .  7  km 
per  hour) ,  but  about  ten  times  as  high  as  Lacombe' s  adjusted  value 
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(0.15  km  per  hour) . 

If  the  net  evaporation  per  square  centimeter  in  the  Mediterranean 
is  about  150  cm/year,  and  all  the  water  deficit  were  to  be  balanced 
by  the  inflow  through  the  Strait  of  Gibraltar,  a  maximum  net  inflow 
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of  0.4  5  km  /hour  would  be  required.  Furthermore,  heat  budget 
computations  at  FNWC  for  the  Mediterranean  require  a  larger  water 
exchange  than  heretofore  estimated.  More  accurate  insight  into 
these  problems  is  expected  from  the  two-layer  HN  model. 

7.  SUMMARY. 

a.  The  HN  model  with  two  open  boundaries  can  be  used  for 
computation  of  tides  and  currents  (including  wind  currents)  in 
relatively  deep  straits. 
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b.  The  tidal  phases  and  amplitudes  computed  with  the  HN 
model  for  the  Strait  of  Gibraltar  verify  well.  The  maximum  differ¬ 
ence  has  been  30  minutes  in  time  of  high  water  and  20%  in  amplitude. 

c.  The  net  (or  rest)  tidal  current  gives  a  net  inflow  into  the 
Mediterranean.  The  inflowing  tidal  currents  are  stronger  and  also 
last  somewhat  longer  than  the  outflowing  tidal  currents. 

d.  The  computed  net  inflow  is  somewhat  larger  than  previously 
available  estimates.  However,  heat  budget  computations  require  a 
large  water  exchange.  Multilayer  HN  models  will  shed  further  light 
on  this  problem. 
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Figure  10  Currents  (tidal  and  wind)  in  Bay  of 

Gadiz  at  low  water  at  Tarifa,  computed 
with  HN  model  (wind  6m, sec  ^  from  E). 


UF  GIBRALTAR  (LOCATIONS  SEE  FIG,  12,  AFTER  LACOHBE  1961), 
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FIGURE  14.  TkAiJSPORT  OF  WATER  THROUGH  STRAIT  OF  GIBRALTAR 
FROM  SURFACE  TO  IOOm.  (HN  MODEL) 


